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基于 改进 Qsplat 算法 的 三 维 动态 积 云 模拟 研究 
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(辽宁 工程 技术 大 学 电子 与 信息 工程 学 院 , 辽宁 芯 芦 岛 125105) 


摘 要: 尽管 人 们 对 云 模拟 进行 了 大 量 研 究 ， 但 是 现 有 的 云 模 拟 算 法 很 难 在 保证 实时 性 的 情况 下 得 到 真实 感 较 强 的 云 
图 形 。 针 对 这 个 问题 ， 提 出 一 种 结合 Qsplat 算法 与 IFS 算法 的 粒子 系统 的 积 云 模拟 方法 ， 并 考虑 了 积 云 的 结构 特点 。 
首先 ， 基 于 Qsplat 算法 的 细节 层次 树 型 结构 建立 树 型 结构 粒子 系统 ， 加 速 粒 子 系统 搜索 速度 ; 然后 ， 对 距离 视点 较 近 
的 粒子 使 用 IFS 迁 代 函数 系统 算法 增加 层次 包围 球 树 的 层次 以 丰富 云 体 局 部 细节 来 建立 云 的 外 形 ; 最 后 ， 使 用 球形 图 
元 映射 二 维 纹理 进行 云 的 绘制 。 最 终 实现 了 三 维 动态 积 云 的 快速 实时 模拟 。 实 验 结果 表明 ， 该 方法 能 够 快速 地 生成 实 
时 三 维 动态 积 去 图形， 图 形 真实 感 强 。 
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Research on three-dimensional dynamic cumulus simulation based on improved Qsplat 


Peng Yanfei, Mi Yang, Zi Lingling, Zhao Quanying 
(School of Electronic & Information Engineering Liaoning Technical University, Huludao Liaoning 125105, China) 


Abstract: Although people have done a lot of research on cloud simulation, the existing cloud simulation algorithm can hardly 
obtain the cloud image with strong realism while ensuring the real-time performance. According to this problem, this paper 
proposes a cumulus-based simulation method for particle systems that combine Qsplat with IFS. First, build a tree-structured 
particle system based on the detail-level tree structure of Qsplat for speeding up particle System search speed. Then, for these 


particles which closer to the viewpoint , used the iterative function system algorithm to increase the level of enclosing the ball 


tree to enrich the local detail of the cloud to establish the shape of the cloud. Finally, Finally, rendered the cloud by using two- 


dimensional texture mapping of spherical entities, which realizes the fast real-time simulation of dynamic cumulus. The 


experimental results show that the method can generate real-time three-dimensional dynamic cumulus images quickly. 
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积 云 模拟 的 难点 所 在 。 在 所 有 现存 的 云 模拟 方法 中 ， 基 本 分 为 
两 类 。 一 种 是 基于 物理 模型 的 模拟 方法 ， 典 型 的 有 ，Harris 等 

不 规则 模糊 物体 的 模拟 长 时 间 以 来 是 计算 机 图 形 学 研究 的 ”通过 在 N-S 方程 中 引入 热力 学 方程 ， 使 用 三 维 平 铺 纹理 技术 实 
主要 热点 和 难点 。 云 是 一 种 广泛 存在 的 自然 不 规则 模糊 物体 ， 现 对 云 的 真实 感 模拟 证 ， 唐 勇 、 赵 丽 文 等 人 则 通过 在 N-S 方程 
无 论 是 影视 或 动漫 制作 还 是 游戏 创作 抑或 是 飞行 场景 仿真 等 等 ”中 添加 粘性 项 与 涡 旋 ， 同 时 考虑 温度 、 密 度 等 物理 特性 ， 结 合 
领域 都 无 法 绕 开 这 个 问题 ， 高 效 而 又 逼真 的 云 模拟 技术 对 图 形 ” 单项 散射 光 模 型 对 三 维 动态 云 进行 模拟 中 ， 何 晓 眶 ， 陈 雷霆 等 
学 发 展 和 现代 影视 游戏 工业 以 及 军工 及 航空 业 都 有 着 重要 的 意 。 人 通过 简化 N-S 方程 建立 云 的 运动 模型 ， 同 时 利用 Impostor 技 
义 。 自 然 界 云 的 种 类 众多 , 但 很 多 都 不 常见 , 在 大 多 数 情况 下 ， 术 配 合 Harris 的 云 光照 数学 模型 加 速 泻 染 大 规模 的 云层 图 形 欠 。 


完全 模拟 所 有 种 类 没有 必要 。 积 云 是 自然 界 最 常见 也 是 特征 最 ”这 些 物理 模型 模拟 方法 往往 能 取得 较 好 的 真实 感 去 图形， 但 是 
明显 的 云 ， 完 成 积 云 的 模拟 基本 上 可 以 满足 绝 大 部 分 需求 。 因为 运用 了 大 量 的 物理 学 原理 ， 导 致 运算 过 程 非常 复杂 ， 运 算 
耻 云 最 明显 的 特征 有 以 下 几 种 : 结构 紧密 ， 外 形 不 规则 ， 量 空前 巨大 ,对 计算 机 硬件 的 要 求 过 高 , 并 不 适合 于 实时 模拟 ， 
“花椰菜 ” 状 ; 云 块 扳 立 ， 具 有 明显 的 自 相 似 性 ， 大 部 分 云 块 上 且 由 于 当前 物理 模型 问题 ， 这 些 模拟 方法 对 云 体 的 运动 控制 


上 
= 
边界 明显 , 存在 一 些 云 块 破裂 , 边界 较为 模糊 中。 这些 特 点 就 是 


蔷 
对 
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FE 问题。 另 一 种 模拟 方法 就 是 基于 个 体 生长 模型 的 云彩 模 
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拟 方法 ， 包 插 很 多 近似 模拟 方法 。 例 如 次 怀 清 、 刘 浩瀚 等 人 对 ”小 于 等 于 阔 值 ， 则 不 再 向 下 遍历 并 绘制 该 粒子 329。 
这 种 方法 进一步 改进 ， 通 过 利用 改进 的 Gardner 椭 球 模型 结合 分 形 几 何 是 近代 数学 发 展 的 一 个 重要 分 支 ， 它 是 研究 非 规 
钱 加 纹理 的 方法 实现 了 实时 的 云 模 拟 罩 。 基 于 粒子 系统 的 模拟 。” 则 几何 形态 的 几何 学 。 不 规则 形态 在 自然 界 广 泛 存在 ， 因 而 分 
方法 一 直 认 为 是 不 规则 物体 模拟 领域 中 最 成 功 的 方法 。 潘 秋 羽 。” ” 形 几 何 同时 被 称 为 是 描述 自然 界 的 几何 学 。 普 遍 认为 1967 年 
等 人 利用 粒子 系统 结合 纹理 映射 的 方法 ， 并 利用 Impostor 等 技 。” Mandelbrot 发 表 的 论文 《英国 海岸 线 有 多 长 》 是 分 形 几何 学 正 
术 加 快 绘制 , 从 而 实现 了 三 维 动态 云 的 实时 模拟 铝 , Elhaddad 利 式 产生 的 标志 。 然 而 事实 上 在 此 之 前 分 形 几 何 的 许多 重要 问题 
用 分 子 动力 学 模拟 模型 中 实现 了 三 维 云 的 模拟 。 这 些 方法 一 般 ”早已 提出 , 例如 斯 宾 塞 三 角 、 康 托 尔 集 等 .迭代 函数 系统 (iterated 
能 取得 较为 通 真 的 效果 和 较 好 的 控制 能 力 ， 但 是 普遍 存在 的 问 ”function system, IFS) 是 分 形 几何 学 中 一 个 重要 且 最 有 发 展 前 景 
题 是 当 模 拟 物体 精度 要 求 较 高 时 往往 需要 大 量 的 粒子 ， 这 对 模 ”的 分 支 。 友 代 函 数 系统 将 待 生成 的 图 形 看 做 是 由 许多 与 整体 相 
拟 的 效率 是 严峻 的 考验 。 还 有 的 就 是 传统 的 基于 动态 网 格 的 模 似 的 〈 自 相似 ) 或 经 过 一 定 变 换 与 整体 相似 的 《〈 自 仿 射 ) 小 块 
拟 方法 ,例如 范 晓 舌 , 张 立 民 等 人 利用 N-S 方程 构建 动态 网 格 ， 拼 贴 而 成 ， 通 过 对 原始 图 形 经 过 缩小 、 平 移 、 旋 转 等 操作 绘制 
结合 多 向 散射 模型 ， 同 时 采用 GPU 加 速 的 云 模拟 方法 图 以 及 出 新 的 图 形 ， 和 迭代 函数 系统 生成 的 图 形 永 远 具有 自 相似 性 。 如 
Yusov 等 提出 的 网 格 模 型 外。 这些 方 法 的 模拟 速度 较 快 , 但 是 建 ” 式 (1) 所 示 就 是 一 个 二 维 的 仿 射 变换 。 
模 过 程 非常 复杂 ， 当 模拟 云 体 精度 较 高 时 ， 建 模 的 难度 更 是 极 了 
网 生生 
C f 


有 具 上 升 ， 同 时 还 存在 的 问题 是 ， 模 拟 的 云 体 是 中 空 的 ， 实 现 动 
态 效果 的 难度 很 大 , 容易 导致 模拟 失真 .除了 两 大 类 研究 方向 ， 2 20 
还 有 一 些 其 他 方法 ,例如 考虑 用 户 的 需求 生成 云 模拟 模型 方法 ， ”其 中 : W 是 一 个 二 维 仿 射 变换 ，a，b，c，d，e，f 为 变换 的 参 
Dobashi 结合 用 户 需 求 和 真实 云层 图 片 生成 云 的 模拟 图 形 04504， ” 数 ， 通 过 更 改 这 些 参数 ， 就 可 以 得 到 不 同 的 子 图， 以 满足 模拟 
Goswami 等 人 引入 用 户 设置 的 拉 格 朗 日 三 维 云 模拟 框架 0; 的 要 求 P0。 
Rui 等 结合 气象 观测 数据 与 简单 物理 模型 的 三 维 云 模拟 方法 [3 。 基于 点 绘制 的 Qsplat 绘制 算法 由 于 预 处 理 阶段 较为 复杂 ， 
且 近 些 年 来 , GPU 的 迅速 发 展 为 动态 云 模拟 提供 了 新 的 研究 。 通常 应 用 于 大 规模 静态 点 云 数据 的 处 理 ， 不 适合 动态 点 绘制 。 
思路 。GPU 高 速 并 行 的 特点 ， 对 物理 模型 计算 和 碰撞 检测 都 有 ”和 迭代 函数 系统 迭代 函数 构造 简单 ， 便 于 计算 机 执行 ， 对 自然 景 
非常 明显 的 加 速效 果 041!17。 例 如 Barbosa 等 人 提出 的 基于 坐标 观 的 模拟 往往 能 取得 很 好 的 效果 。 通 过 结合 两 种 算法 ， 可 以 对 
的 流体 模拟 技术 结合 GPU 加 速 计算 技术 方法 [51。 传统 Qsplat 算法 进行 改进 的 到 融合 IFS 的 动态 的 Qsplat 算法 对 
针对 积 云 的 特点 要 求 以 及 现 有 算法 存在 的 不 足 ， 本 文 提 出 “动态 点 云 数 据 进 行 处 理 。 
了 一 种 0 a 2 ”三 维 动态 积 云 模拟 
Qsplat 算法 层次 包围 球 树 结构 粒子 系统 ， 结 合 和 迭代 函数 系统 算 
法 对 层次 包围 球 树 树 的 局 部 粒子 进行 分 形 处 理 ， 得 到 细节 层次 针对 积 云 结构 紧密 、 孤 立成 块 以 及 自 相 似 性 等 特点 ， 本 文 
更 高 的 粒子 数据 ， 最 终 利用 球形 图 元 映射 二 维 纹理 绘制 出 积 云 。” 提出 了 基于 Qsplat 的 云 模拟 算法 。 该 算法 使 用 基于 Qsplat 算法 
图 形 。 经 过 实验 ， 算 法 在 真实 感 和 绘制 速度 方面 取得 了 较 好 的 。” ”的 特殊 树 结构 粒子 系统 构建 云 体 模型 ， 然 后 利用 IFS 迭代 函数 
平衡 。 系统 丰富 云 体 局 部 细节 ， 最 后 利用 球形 图 元 映射 二 维 纹理 对 粒 
子 进行 绘制 。 算 法 的 流程 如 图 1。 
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1 ”Qsplat 方 法 及 IFS 算法 的 原理 


初始 化 树 形 粒子 系统 结构 


2000 年 , 加利福尼亚 大 学 图 形 学 实验 宇 提 出 了 一 种 改进 的 根据 森 形 科 子 系统 结构 建立 基础 Qsp at 层 次 包围 球 树 
大 规模 点 云 数据 加 速 处 理 算法 一 一 Qsplat 算法 。 这 是 一 种 基于 计 看 税 子 变 到 帮 用 三 
点 的 绘制 算法 ， 相 比 于 splat 等 传统 基于 点 的 算法 , 该 算法 具有 ee 
更 快 的 绘制 速度 。 该 算法 的 原理 如 下 。 Re 
Qsplat 的 预 处 理 ， 通 过 长 轴 二 分 法 等 方法 对 点 云 数据 进行 人 
划分 处 理 ， 将 点 云 树 构建 为 层次 包围 球 的 树 状 结构 。 该 树 结构 和 
包含 了 不 同 细节 层次 的 点 模型 ， 树 的 根 节点 代表 整个 模型 ， 树 
中 的 各 层次 节点 代表 不 同 精度 下 的 点 模型 采样 ， 每 个 叶子 节点 a 
代表 一 个 或 多 个 点 云 模型 中 的 点 数据 。 树 中 各 节点 最 为 一 个 基 ee 
图 1 本 文 算法 流程 图 
础 绘制 图 元 ， 其 属性 包括 位 置 、 半 径 、 法 线 以 及 颜色 等 。 树 中 
节点 的 属性 计算 要 自 底 向 上 进行 。 在 包围 球 树 结构 建立 后 ， 进 ” 2.1 三 维 动态 积 云 的 建 模 
行 遍历 绘制 。 人 遍历 标准 :从 根 节点 开始 ， 判 断 当前 节点 包围 球 2.1.1 积 云 粒子 系统 的 树 型 结构 
的 屏幕 投影 是 否 大 于 阔 值 。 若 大 于 阔 值 ， 则 继续 向 下 遍历 ， 若 积 云 的 结构 紧密 、 孤 立成 块 且 运动 复杂 。 粒 子 系统 方法 是 
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最 适合 模拟 不 规则 模糊 动态 物体 的 方法 ， 但 是 简单 的 粒子 系统 
在 控制 方面 存在 问题 。 针 对 这 些 问 题 ， 为 了 得 到 真实 感 强 且 易 
于 控制 的 积 云 模型 ， 并 方便 构建 动态 的 Qsplat 细节 层次 树 ， 本 
文 提 出 一 种 特殊 的 树 型 组 织 结构 来 构建 粒子 系统 ， 并 在 具有 父 
子 关 系 的 粒子 间 添 加 弹力 来 约束 粒子 的 运动 范围 ,如 图 2 所 示 。 


“1 


图 2 ” 树 型 结构 粒子 系统 示意 图 


算法 首先 确定 根 节 点 粒子 的 位 置 与 半径 (也 就 是 云 体 的 位 
与 大 小 )， 然 后 在 根 粒子 的 半径 内 创建 多 个 子 节点 粒子 ， 每 一 
个 子 节 点 再 创建 自己 的 七 又 树 子 树 ， 多 个 七 又 树 通 过 最 终 连接 
成 一 个 整体 。 整 个 粒子 树 中 根 节点 只 有 在 需要 移动 整个 云 体 是 
会 产生 位 移 ， 除 此 外 不 受 任何 力 影 响 ， 其 位 置 固定 ， 不 会 随时 
间 变 化 而 移动 。 粒 子 系统 中 除根 节点 外 的 粒子 都 会 受到 重力 、 
风力 、 空 气 阻力 的 影响 。 其 中 值得 注意 的 是 ， 除 了 根 粒子 外 的 
粒子 所 受 弹 力 仅 来 自 于 它 的 父 节点 粒子 ， 这 样 可 以 更 好 地 控 和 
整个 粒子 系统 的 体积 ， 粒 子 系统 中 的 粒子 《除根 粒子 ) 都 有 自 
己 的 半径 ， 半 径 值 会 随时 间 和 粒子 位 置 的 变化 而 发 生 改变 。 
2.1.2 积 云 粒 子 的 运动 方程 

队 云 的 运动 主要 表现 为 由 最 初 的 小 淡 积 云 逐渐 向 上 、 向 四 
司 生 长 ， 形 成 各 种 结构 紧密 的 不 规则 的 花椰菜 形 。 基 于 物理 方 
法 的 模拟 往往 具有 较 大 的 运算 量 ， 对 计算 机 硬件 具有 很 高 的 要 
求 ， 不 利于 大 规模 实时 模拟 。 因 此 ， 本 文 对 粒子 运动 模型 进行 
了 简化 。 

根 节点 粒子 作为 销 点 确定 整个 云 体 的 位 置 与 大 概 范围 ， 各 
层次 的 子 节点 粒子 收 风 力作 用 产生 运动 ,模拟 积 云 生长 的 过 程 ， 
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模拟 积 云 向 四 周 扩张 以 及 向 上 生长 的 特点 。 模 型 
函数 表达 式 如 式 (2)。 
/ixz /ixx 
yn 一 XZ 
”六 J 


其 中 ,入 为 风力 在 x-z 平 面 方 向 分 量 的 模 ,表达 式 如 式 (3) 所 示 。 


的 风 场 模型 ， 
的 风力 是 积 云 生长 的 动力 ， 


XXX 十 万 xy 十 


矿 =6.25x104x2+z2+0001xy+0.1 G) 
包 为 风力 在 y 方 向 分 量 的 模 ， 表 达 式 如 式 (4)。 


户 =-1.43x10” Vx +z -1.43x10™* 7+0.004 


x、y 和 z 为 三 维 空间 坐标 ， 通 过 控制 常数 系数 的 值 可 以 改 
变 积 云 在 水 平方 向 与 垂直 方向 生长 的 速度 。 

粒子 所 受 空气 阻力 为 其 运动 的 阻力 之 一 ， 
(5) 所 示 。 


函数 表达 式 如 式 


F,,, =—1.5(speedx + speedy + speedz) (3) 


其 中 ，-1.5 为 空气 阻力 系数 值 ， 空 气 阻力 系数 值 越 大 ， 粒 子 速 
度 下 降 越 快 。speedx、speedy 和 speedz 分 别 为 坐标 x、y 和 z 方 
向 的 速度 分 量 。 

粒子 弹力 是 约束 云 体 的 主要 作用 力 ， 函 数 表 达 式 如 式 (6) 所 


已 =0.05(1.0-d) (0) 


spring 


改变 弹力 系数 的 值 可 调节 积 云 的 
间 的 距离 ， 表 达 式 如 式 


其 中 ，0.05 为 弹力 系数 的 值 ， 
紧密 程度 ，d 为 当前 粒子 与 其 父 粒子 2 


(7) 所 示 。 
全 > 2 
d=A\J(z-z) +(y—y) +(z—z,) (7) 
x、y、 以 及 z 为 当前 粒子 位 置 坐标 ，x0、y0 和 z0 是 其 父 


节点 粒子 位 置 坐标 。 
2.1.3 结合 IFS 算法 的 改进 Qsplat 算法 
单纯 的 的 粒子 系统 为 达到 较 高 的 真实 感 往 往 会 需要 大 量 的 


在 粒子 与 其 父 节 点 粒子 间 添 加 弹力 约束 粒子 运动 范围 、 控 制 云 
体 的 生长 不 至 于 脱离 范围 导致 云 体 空洞 。 本 文 模型 中 除根 节点 
粒子 外 所 有 粒子 运动 过 程 中 受 重力 、 来 自 父 粒 子 的 弹力 、 风 力 
以 及 空气 阻力 的 作用 ， 县 体 情 况 如 图 3 所 示 。 其 中 ，G 是 粒子 
所 受 重 力 ，spring 表示 弹力 ，F 表示 风力 ，f 表 示 空 气 阻力 。 


a 


F 
Lo 
6 Z 
Spring, 


图 3 积 云 粒子 受 力 分 析 图 


> 


为 了 高 效 真 实地 模拟 积 云 的 生长 效果 ， 本 文 提出 一 种 简化 


粒子 来 提高 模拟 图 形 的 精度 ， 然 而 这 又 会 降低 模拟 的 效率 。 为 
了 解决 这 个 矛盾 ， 本 文 提 出 了 一 种 改进 的 Qsplat 算法 来 保证 在 
不 降低 模拟 图 形 的 精度 的 条 件 下 的 模拟 效率 。 

Qsplat 算法 一 般 被 用 来 加 快 处 理 大 规模 的 静态 点 云 数 据 ， 
它 需 要 经 过 复杂 的 预 处 理 阶 段 ， 这 不 利于 动态 实时 图 形 绘制 。 
因此 本 文 对 Qsplat 算法 进行 改进 ， 固 化 具有 一 定 层 数 的 基础 细 
节 层 次 树 简 化 预 处 理 ， 来 适合 本 文 算法 的 动态 实时 的 需要 。 同 
时 ， 算 法 为 了 提升 模拟 图 形 局 部 细节 层次 ， 需 要 为 一 些 距 离 视 
点 较 近 的 粒子 扩展 细节 层次 ， 因 此 在 基本 Qsplat 算法 基础 上 融 
合 IFS 和 迭代 函数 系统 ， 对 这 些 粒子 进行 分 形 途 代 处 理 ， 增 加 局 
部 粒子 ， 丰 富 模拟 图 形 局 部 层次 细节 ， 提 高 图 形 的 真实 感 。 

对 于 给 定 的 需要 分 形 处 理 的 叶子 节点 M(x,y,z)， 
线性 变换 产生 i 个 点 Mi(xi,yi,zi),i=0,1,…,7。 

线性 变换 的 表达 式 如 式 (8)。 


DH 


经 过 一 次 
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所 示 的 纹理 图 片 来 模拟 积 云 外 形 。 


X, a b cllx m 
de flyl+t+in (8) 


ll 


2; sg h jillz 1 


为 使 迭代 收 僵 ， 每 一 次 迭代 必须 是 压缩 变换 。 设 根 节 点 半 
径 为 R， 当 前 待 分 形 处 理 粒子 半径 为 r， 压 缩 因子 C 的 表达 式 


图 6 积 云 粒子 纹理 


如 式 (9)。 
C=R/r (9) L 体 的 绘制 步 又 如 下 : 

联 立 得 如 式 (8) 的 方程 组 ， 求 出 线性 变换 的 相关 系数 ， 就 可 a) 对 所 有 待 绘制 粒子 根据 深度 由 大 到 小 排序 并 存储 至 待 
以 得 到 相应 的 IFS 码 。 绘制 缓存 区 ; 

改进 Qsplat 算法 具体 步骤 如 下 : b) 饥 历 待 绘制 粒子 缓冲 区 ， 为 加 快 绘制 速度 ， 本 文 实验 采 

a) 基 于 树 型 粒子 结构 建立 如 图 4 所 示 具 有 固定 深度 基础 ”用 directx9 内 置 球体 绘制 函数 来 预先 计算 球形 图 元 数据 ， 读 得 
Qsplat 细节 层次 树 ; 粒子 数据 后 直接 在 相应 位 置 绘制 出 球体 ; 

b) 遍 历 Qsplat 细节 层次 树 。 若 当前 粒子 的 半径 小 于 浆 值 ， c) 利 用 纹理 混合 Blend 等 技术 完成 积 云 粒子 绘制 。 


停止 遍历 ， 将 当前 粒子 放 入 待 绘 制 粒 子 缓存 区 ， 否 则 继续 向 下 3 ”实验 与 分 析 
遍历 ， 若 粒子 深度 被 遮挡 且 深 度 超 过 阔 值 ， 停 止 遍 历 ， 若 粒子 
闵 值 大 于 阔 值 但 是 当前 粒子 是 叶子 节点 则 转 到 第 c) 步 ; 本 文 使 用 VS2010 与 Directx9c 实现 本 文 算法 , 所 有 实验 在 
c) 基础 细节 层次 树 的 细 化 分 形 处 理 。 若 当前 粒子 是 叶子 节 ”CPU: Inteli3-2330， 主 频 : 2.4 GHz， 内 存 : 4GB， 显 卡 : NVS 
点 粒子 且 半 径 大 于 阔 值 ， 对 当前 粒子 用 从 代 函 数 系统 进行 分 形 3100M， 分 辨 率 : 1280*1024 的 机 器 上 进行 
处 理 并 生成 它 的 临时 子 树 ， 继 续 检 查 叶 子 节点 ， 若 满足 条 件 则 经 过 实验 ， 发 现 模 拟 最 低 要 求 粒子 树 有 两 层 8 个 粒子 ， 
停止 ， 否 则 继续 执行 c) 步 ; 7 中 所 示 为 模拟 积 云 粒子 数 为 8、Qsplat 树 基础 层次 为 2、 云 体 
b) 步 。 初始 大 小 4km 不 同时 刻 的 云 模拟 图 形 。 
当 粒 子 树 为 三 层 57 个 粒子 时 实验 在 粒子 系统 初始 化 、 模 
加 拟 效果 以 及 运行 帧 率 之 间 取 得 一 个 较 好 的 平衡 。 图 8 中 所 示 为 
模拟 积 云 粒子 数 为 57 、Qsplat 树 基础 层次 为 3、 云 体 初始 大 小 
4km 不 同时 刻 的 效果 图 。 


双 


d) 计算 粒子 的 运动 速度 ， 更 新 粒子 属性 


主 
讽 
五 


图 4 Qsplat 层次 包围 球 结构 树 


2.2 三维 动 态 积 云 的 绘制 

队 云 主要 由 浓密 的 球状 云 泡 组 成 ， 大 量 球形 气泡 堆积 形成 
各 种 不 规则 外 形 。 考 虑 这 个 问题 ， 本 文 使 用 的 粒子 以 球形 图 元 
作为 基础 绘制 单元 。 积 云 粒子 为 球形 图 元 映射 二 维 纹理 图 片 绘 
制 而 成 。 同 时 ， a 避免 图 元 之 间 空 洞 的 出 现 ， 球 形 图 元 的 半 
径 会 随 着 粒子 的 运动 而 变化 ,算法 使 用 的 球形 图 元 如 图 5 所 示 。 


(b)T=10s 


pi 


和 (a)T=15s (b)T=20s 
| ~- 一 一 图 7 粒子 数 为 8 时 积 云 模拟 图 形 
帧 率 是 实时 算法 最 重要 的 指标 之 一 ， 图 9 所 示 为 积 云 粒子 


数 为 8 和 积 云 粒子 数 为 57 时 模拟 帧 率 随时 间 的 变化 曲线 。 
实验 可 知 ， 算 法 运行 的 帧 率 一 直 保持 较 高 的 值 ， 完 全 满足 实时 
粒子 系统 中 的 粒子 分 类 以 距离 视点 的 大 小 为 依据 ， 处 于 视 ”模拟 的 需要 ， 并 且 随 着 时 间 增 加 ， 实 验 的 帧 率 会 逐渐 稳定 。 


图 5 球形 粒子 图 元 示意 图 


点 周围 固定 距离 矩形 区 域内 的 粒子 利用 IFS 算法 进行 进一步 处 结合 效果 图 与 帧 率 -时 间 折 线 图 分 析 可 知 , 本 文 算法 生成 的 
理 ， 最 终 采 用 多 个 球形 粒子 映射 二 维 纹理 图 片 绘制 ， 处 于 该 区 云 生 成 的 云 结构 紧密 ， 较 好 地 模拟 了 积 云 的 自 相似 性 ， 实 现 了 
域外 的 粒子 直接 采用 球形 图 元 映射 二 维 纹理 进行 绘制 。 积 云 的 向 周围 扩张 、 向 上 生长 的 运动 效果 ， 实 现 了 三 维 积 云 的 


考虑 积 云 紧密 的 结构 及 凹凸 不 平 的 外 表 ， 本 文 使 用 如 图 6 ”动态 模拟 。 算 法 在 模拟 粒子 树 只 有 2 层 8 个 粒子 时 模拟 的 积 云 


201804.02023v1 


chinaXiv 


期 
上 二 


ChinaXiy 会 作 期 刊 
录用 稿 彭 曼 飞 ， 等 : 二 下 过 Qo 浊 骂 积 云 模拟 研究 


云 体 粗糙 ， 云 体 紧密 程度 较 差 ， 模 拟 的 平均 帧 率 为 607， 在 运 ,| 

行 超过 30s 后 的 稳定 平均 帧 率 为 311 ; 在 模拟 粒子 树 达到 3 层 [mh 这 粒 于 系统 
57 个 时 ， 模 拟 图 形 云 体 更 加 紧密 ， 云 体 精细 度 较 好 ， 模 拟 的 平 me | 
均 帧 率 为 220， 算 法 运行 30s 后 的 平均 帧 率 为 102。 由 实验 本 | 

知 ， 算 法 可 以 实现 较为 逼真 的 动态 积 云 模拟 ， 并 且 可 以 通过 增 

加 粒子 数量 增加 模拟 云 的 规模 和 细节 质量 ， 当 模拟 粒子 数量 为 上 | 

57 时 ， 算 法 生成 积 云图 形 的 效果 较 好 ， 完 全 满足 实时 要 求 。 nn 


图 10 帧 率 随 粒子 数 变 化 折线 图 
, 本 文 算法 还 与 文献 [5,6] 提 出 的 算法 进行 了 比较 ， 如 表 1 所 
示 是 三 种 算法 的 基本 属性 情况 。 ee 
2 , Gradner 椭 球 模型 结合 纹理 映射 来 实现 云 的 模拟 ， 但 是 没有 实 


现 动 态 效 果 ; 文献 [6] 通 过 粒子 系统 结合 Impostor ee 了 
三 维 动态 云 的 模拟 。 


表 1 本 文 算法 与 文献 [5] 和 文献 [6] 的 对 比 
算法 提出 建 模 方法 三 维 体 云 是 否 动 态 是 否 实时 
文献 [6] ”Gradner 椭 球 模型 是 将 是 
一 文献 [8] 粒子 系统 是 是 


(a)T=5s (b)T=10s 


是 
(c)T= 15s (d)T=20s 本 文 算法 , 型 结 构 是 是 是 
图 8 粒子 数 为 57 时 积 云 模拟 图 形 粒 了 系统 
| -二 二 本 文 算法 与 文献 [5]、 文 献 [6] 都 为 实时 算法 , 为 了 进行 算法 
瑟 粒 子 数 为 8 
joo | |-e 粒子 数 为 57 的 对 比 ， 本 文 对 粒子 数 分 别 为 8、57 和 400 时 三 种 算法 分 别 进 


行 了 实验 ， 运 行 效果 图 如 图 11~13 所 示 (本 文 算法 效果 图 均 选 
取 自 帧 率 稳定 时 实验 效果 图 )。 结 合 图 12、13 可 以 看 出 ， 与 文 
献 [5] 相 比 ， 粒 子 数 分 别 为 8、57 和 400 时 ， 改 进 Qsplat 算法 的 
加 速 比 率 分 别 为 116.55%、380.0%、1850.0%， 同 时 改进 Qsplat 
算法 生成 的 动态 云 拥 有 更 高 的 真实 感 效 果 。 从 实验 结果 可 以 看 
出 本 文 算法 只 需要 少量 粒子 就 可 以 模拟 云 的 动态 效果 ， 相 比 动 
态 云 模拟 的 文献 [3] 算 法 ， 本 文 算法 的 帧 率 有 较 大 提升 ， 实 时 性 

为 验证 改进 Qsplat 算法 的 加 速效 果 ， 将 单纯 粒子 系统 与 结 ” 更 好 ; 结合 图 11 与 图 13 可 以 看 出 ， 与 文献 [6] 的 算法 相 比 ， 虽 
合 改进 Qsplat 算法 的 粒子 系统 的 运行 效率 进行 比较 。 实 验 结果 然 运 行 帧 率 不 如 ， 但 是 仍 能 完全 满足 实时 性 要 求 ， 并 且 相 比 文 
的 帧 率 随 粒子 数 增加 的 折 现 图 如 图 10。 单 纯 粒 子 系统 的 平均 帧 。 献 [6] 算 法, 改进 Qaplat 算法 不 仅 实现 了 动态 积 云 效果 , 而 且 生 
率 为 126.2; 结合 改进 Qsplat 算法 的 粒子 系统 的 平均 帧 率 为 ”成 的 积 云图 形 比 文献 [6] 生 成 的 云彩 图 形 更 具 立 体感 ， 真 实感 更 
802.7。 结 果 表 明 改 进 Qsplat 算法 有 明显 的 加 速效 果 。 强 。 


时 间 /s 
图 9 帧 率 随 时 间 变 化 折线 图 
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(a) 粒 子 数 8， 平 均 帧 率 802 (b) 粒 子 57， 平 均 帧 率 241 (c) 粒子 400， 平 均 帧 率 78 
图 11 文献 [6] 效 果 图 及 相关 信息 和 帧 率 
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FPS:143.719 


(a) 粒子 数 8， 平 均 帧 率 145 


(b) 粒子 数 57， 平 均 帧 率 20 


FPS :18.981 


(c) 粒子 数 400， 平 均 帧 率 4 


图 12 文献 [8] 效 果 


图 及 相关 信息 和 帧 率 


(a) 粒子 数 8, 帧 率 314 


(b) 粒 子 数 57， 帧 率 96 


(c) 粒子 数 400， 帧 率 40 


4 ”结束 语 


职 云 的 不 定型 外 形 给 模拟 带 来 了 很 大 的 阻碍 ， 如 何在 实时 
性 与 模拟 的 真实 感 之 间 找 到 平衡 是 积 云 实时 模拟 的 难点 。 本 文 
提出 了 一 种 基于 动态 Qsplat 算法 的 树 型 粒子 系统 积 云 模拟 算法 
在 保证 模拟 精度 的 条 件 下 尽 可 能 地 提高 模拟 效率 ， 算 法 提出 的 


并 建立 动态 层次 包围 球 树 加 快 绘制 3 丰富 积 云 局 部 细节 。 
实验 结果 表明 ， 算 法 在 保证 模拟 图 形 真实 感 的 基础 上 ， 实 现 了 
三 维 动态 积 云 的 快速 实时 模拟 ， 可 应 用 于 飞行 场景 仿真 、 游 戏 
制作 等 领域 。 

本 文 没有 对 云 的 光照 问题 进行 深入 研究 ， 因 此 无 法 实现 真 
实 光照 条 件 下 云彩 色彩 的 变化 ， 改 进 Qsplat 算法 的 结构 特点 限 
制 了 该 算法 的 运用 场景 ,只 能 用 于 云 等 运动 不 剧烈 的 物体 模拟 。 
这 些 问题 都 是 后 续 研 究 的 重点 
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